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Abstract. We propose a general method for deriving asymptotics of the Maclaurin 
series coefficients of algebraic functions that is based on a procedure of K. V. Safonov 
and multivariate singularity analysis. We test the feasibility of this this approach by 
experimenting on several examples. 



I. Introduction 

Let N + denote the set of (strictly) positive natural numbers. For m G N + , f3 G N m , 
and x G C m , let Xi denote component i of x, let (3\ = (3\\ • • ■ (3 m \, n{3 = (n/?i, . . . , n/3 m ), 
x 13 = xf 1 ■ ■ ■ x 1 ^ 1 , and d 13 = <9f 1 ■ ■ • df^ 1 , where dj is partial differentiation with respect to 
coordinate j. For a domain H C C m (an open connected set) containing the origin and a 
holomorphic function / : H — >■ C, let fp — d^f(0)//3\ for G N m , the coefficients of the 
Maclaurin series for /. 

In |SafOO] K. V. Safonov proved the following theorem. 

Theorem 1.1. If / :C C d_1 — > C is an algebraic function holomorphic in a neighborhood 
of the origin and every term of the Maclaurin series of f(x) contains a factor of x\, then 
there exists a rational function F :C C d — > C holomorphic in a neighborhood of the origin 
and a matrix M G GLd_i(Z) with nonnegative entries such that the following relation 
between the Maclaurin series coefficients of / and F holds: 

fp = ^7,71 > 

where (3 G N 6 ^ 1 and 7 = f3M. 

Loosely speaking, every (d — l)-variate algebraic power series is a generalized diagonal 
of a ci-variate rational power series. Moreover, Safonov gave a constructive procedure for 
finding F given an annihilating polynomial for /. 

In this report we propose a new approach to deriving asymptotics of Maclaurin series 
coefficients of algebraic functions /, namely, (1) apply Safonov's procedure to / to get 
a multivariate rational function F, then (2) apply multivariate singularity analysis to 
F to derive the asymptotics of the Maclaurin series coefficients of F in the direction 
appropriate to /. In this way we reduce deriving asymptotics of algebraic functions to 
deriving asymptotics of rational functions. This is promising, because current general 
methods for the former task exist only in the univariate case (see |FS09[ Chapter VII] 
for a summary of these), whereas current general methods for the latter exist in the 
multivariate well (see |PW08] for a summary of these). 

We test the feasibility of this approach by experimenting on several specific algebraic 
functions. In some cases we meet with difficulty because of gaps in the present theory of 
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multivariate singularity analysis for rational functions. We note these gaps and indicate 
directions for further research. 



2. Examples 

Let us apply our new approach to several specific algebraic functions, beginning with 
univariate ones. Their asymptotics are already known, so we can check our answers. We 



used Maple 11 to do the calculations, and our worksheets are available at http : / / www . cs . auckland . ac . 

A description of Safonov's procedure, which employs resolution of singularities, can 
be found in |SafO O] . and a description of multivariate singularity analysis for rational 
functions can be found in |RW12j and [PW08J. We will assume you are familiar with this 
material. Warning: |SafOO] contains some confusing passages; see Section [3] for a helpful 
annotation. 

Example 2.1 (Catalan numbers right-shifted). Success. 
Consider the algebraic function 



i-yr^4^ 

f(xi) = — 



2 

It is the generating function of the Catalan numbers multiplied by scj, and so its sequence 
of Maclaurin coefficients is the Catalan sequence prefixed by 0. We seek asymptotics for 

fn- 

Since the generating function of the Catalan numbers is holomorphic in a neighborhood 
of the origin, we may apply Safonov's procedure to /. Doing so, we get that f n = F n>n , 
where 

Xi + x 2 — 1 

So the Maclaurin coefficients of / are the diagonal Maclaurin coefficients of F. Notice 
that this is always the case for univariate / since d = 1 implies M — (1) in Theorem ll.il 
Notice also that in this example / is combinatorial, that is, its Maclaurin coefficients 
are all nonnegative, but F is not combinatorial. 

To derive asymptotics for f n , we derive asymptotics for F na with a = (1, 1). To do this 
we apply multivariate singularity analysis to F. Letting H equal the denominator of F, 
we find that V = {x G C 2 : H(x) = 0} is globally smooth. There is one critical point of 
F for a, namely c = (1/2, 1/2), and c is clearly strictly minimal. Thus by the formulas 
of |RW12j ( |RW12| formulas), 



F = 4 n 

± na ^ 

as n — > oo, which is correct. 



_40Fn 3 / 2 32^n 5 / 2 



+ 0{n- 7 ' 2 ) 



Example 2.2 (Supertrees). Failure. Critical points not minimal. 
Consider the polynomial 

P(y, x x ) = y A - 2y 3 + (1 + 2x 1 )y 2 - 2y Xl + Ax\. 

It is an annihilating polynomial for the combinatorial generating function f{x\) for 
supertrees (see |FS09] for more details). The first few Maclaurin coefficients of / are 
known to be 0, 0, 2, 2, 8, 18, 64, 188. Applying Safonov's procedure to / via P we get that 
fn = F n,n, where 

, 2xiXn(2xnX 2 — 3xnXi + x 2 + 2x^xx — 1) 

F(Xi,X2) = 3 9 • 

x%x{ — 2x%X\ + x 2 + 2x 2 X\ — 2 + 4xi 
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To derive asymptotics for f n , we derive asymptotics of F na where a = (1, 1). Letting 
H equal the denominator of F, we find that V is globally smooth. There are three critical 
points of F for a, namely (3/16 - (1/16)^, 1 + v^), (3/16 + (1/16)^, 1 - v 7 ^), and 
(1/8, 2). However, none of these points is minimal. To see this, for each of these points q 
we compute the first nonzero term of the Maclaurin series of g (which depends on q; see 
[RW12j for a definition of h and g) and get at 4 with a < 0. Thus 3tg < near 0, hence 
\h(t)\ < \q 2 \ near 0, hence there exists a point q' G V such that \q[\ = |gx| and \q' 2 \ < \q' 2 \, 
hence q is not minimal. So the [RW12j formulas are not guaranteed to work. 

Plugging c = (1/8, 2) into the formulas anyway we get 



4" 



' ' +0(n- 7 / 4 ) 



_8r(3/4)n 5 /4 

which is nonreal and thus incorrect since f n >0. The correct answer is 

1 



4" 



_8r(3/4)n 5 /4 
(sec [FS09, Chapter VI] for instance). 



+ O (n- 7 / 4 ) 



Example 2.3 (Basic univariate). Failure. No critical points. 
Consider the simple combinatorial algebraic function 

ft \ Xl 



Applying Safonov's procedure to / we get that /„ = F n>n , where 

x 2 Xi(—3x 2 + 3x%Xi — 2 + 2x 2 Xi + X\ — 2x\ + 2x2^1) 



F(x 1 ,x 2 ) 



-x 2 + x\x\ — 2 + 2x 2 X\ + X\ 



To derive asymptotics for f n we derive asymptotics for F na where a = (1, 1). Letting 
H equal the denominator of F, we find that V is globally smooth. However, there are no 
critical points of F for a. Thus we can not apply the |RW12] formulas. 

We try to get around this as follows. When a = (1, A) for < A < 1, there are two 

critical points of F, namely, (|A 2 + (|A - |)\/A 2 — 10A + 9-A + f , 3 ~ 3A+ 2 ff_y° A+ g ) and 

(|A 2 - (f A - §)VA 2 -10A + 9 - A + f, 3 - 3A -ff- 10A+ -)- Plugging the second and first 
critical points into the |RW12j formulas and taking the limit as A — > 1~ we get 

Fna / ~\~ OO, 

V27m 

and its negative, respectively. So both expansions are clearly incorrect, but the leading 
term of the first is curiously close to that of the correct expansion, which is 

F na = -j= + (n- 3 / 2 ) 



(see [FS09j Chapter VI] for instance). 

Alternatively, searching for critical points of F for a = (1,1) in CP 2 yields two points, 
namely (0,1,0) and (0,0,1) in homogeneous coordinates. However, mapping these back 
down to C 2 eventually leads to a trivariate generating function with cone point singulari- 
ties, a type of singularity for which there is presently no multivariate singularity analysis 
worked out. 

3 



Example 2.4 (Basic bivariate). Success. 
Consider the bivariate algebraic function 



f(x 1 ,X 2 ) = X 1 y/l -Xi- X 2 , 



which appears in |SafOO|. Example 2]. Following Safonov's procedure we get that f n p 
F n f3M where M — ( \ \ ) and 



F(x 1 ,x 2 ,x 3 ) 



x 3 xi(3x 3 + 2 + xi + x 2 X\ + 2x\ 
x 3 + 2 + xi + x 2 xi 



Letting H equal the denominator of F, we find that V is globally smooth. There is one 



critical point of F for 0M, namely 



-/3i p 2 



■1). Since F is not combinatorial, we have 



no shortcuts for checking whether or not this point is minimal. Proceeding anyway and 
plugging the critical point into the [RW12J formulas we get 



n/3M 



01 



02 



01 



7T(A + 2 ) 2 VM~2U 2 



+ O (n~ 3 ) 



as n — > oo uniformly as 0M normalized to length 1 varies within a compact set in 
This expansion agrees with experimental checks. 

Example 2.5 (Another basic bivariate). Failure. No critical points. 
Consider the bivariate algebraic function 



f(x u x 2 ) 



X\ 



y/l — X 1 — X 2 

Following Safonov's procedure we get that f n p = F n pu where M = ( { ° J and 



F(xi,X2,X3) 



^3^l( — 3^3 + 3x^X1 + 3x^xiX2 — 2 + 2x3x1 + 2x3x12:2 + xi + X1X2 — 2x| + 2xijXi + 2x|xiX2) 
— X3 + x|xi + X3X1X2 — 2 + 2x3x1 + 2x3x1x2 + xi + X1X2 



Letting H equal the denominator of F, we find that V is globally smooth. However 
there are no critical points of F for 0M. Thus it is unclear how best to proceed in deriving 
asymptotics. It seems we should first understand the no-critical-points phenomenon in 
the simpler bivariate setting of Example 12.31 

3. Comments and Questions 

As we see, to implement our new approach for deriving asymptotics of Maclaurin coef- 
ficients of algebraic functions in full generality, several gaps in the theory of multivariate 
singularity analysis of rational functions need to be filled, namely, 

• How does one check for minimality in the non-combinatorial case? 

• How does one handle non-minimal critical points? 

• How does one handle the case of no critical points? 

Thus more research is required. 

We end with an observation on Safonov's procedure and a short annotation of |SafOO] . 

Proposition 3.1. Let F be the rational function from Theorem ll.il At a critical point 
of F for a := (mi, . . . , rrid-i, mi) G whose last component is nonzero, the numerator 
of F vanishes. 
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Proof. From Safonov's procedure, F(x) has the form "polynomial in x + (power of x) ■ 
R(x)" , where 

R{x,y) = -, — ^7 x— , 

k P(xiy,x 2 ,...,x d - 1 ,y) 

y = Xd (for easy reading), k e N+, P is a polynomial, and the variable substitution is 
understood to occur after the partial differentiation. By definition, the denominator of 
a rational function vanishes at a critical point. Thus by the specific form of F above, it 
suffices to prove the proposition for the rational function R(x,y). 

Let J be the denominator of R. Among the critical equations of R for a are 

J (x, y) = 
a^xxdiJix.y) = a^ydaJ&y), 

that is, 

P(x 1 y,x 2 , . . .,x d ^x,y) = 
xxy(diP)(xxy,x 2 , ■ ■ . ,x d -!,y) = y [x 1 (diP)(xiy,x 2 , ■ ■ .,x d -!,y) + {d d P)(x 1 y,x 2 , ■ ■ ■ ,x d -i,y)} . 

Thus at a critical point where y ^ 0, P(xiy, x 2 , . . . , x d -i, y) — and (d d P){x\y, x 2 , ... , x d -i, y) 
0, as desired. □ 

Corollary 3.2. For a critical point of F whose last component is nonzero, the first term 
in the [RW12j formulas for the asymptotic expansion of F is zero. 

Proof. By [RW12|. Theorem 3.4], the first term in the |RW12j formulas for the asymptotic 
expansion of F na contains a factor of the numerator of F evaluated at the critical point. 
Thus by the previous proposition this term is zero. □ 

In retrospect we should have expected no critical points in Example I2.3[ for in the 
presence of a critical point the [RW12J formulas and Corollary 13.21 predict no rC 1 ! 2 term 
in the expansion for F na whereas the correct expansion does have an rT 1 ! 2 term. 

Annotation to [SafOOj 

• page 261, line 4 of the main body of the article: z should be z\ 

• page 263, lines 2-5: it appears that Safonov is defining his term 'A-diagonal' here 

• 264, 6: R should be R(w, z) 

• 266, 1: 'of the degree g' should be 'of some degree q' 

• 266, 7: ( mn should be £™ n and the remark 'for some mi, . . . ,m n e N' should be 
added 

• 267, 12: the asterisks should be to the left of the rightmost parenthesis 

• 269, 9: w(C) should be u(v,() 

• 269, last line: here Safonov states that there exists such an I but does not explain 
how to find it; it seems one need to look at a few terms of the Maclaurin expansions 
of both h and a and find the I that works by trial and error; note that I is needed 
for the definition of R^ on page 270 

• 270, 2: Ci nl should be (" ln 

• 270, 11: fri(C) should be /i (1) (C) 

• 270, 13: the leftmost occurrence of R^'(z ,z) should be R^ 2 '(z ,z) 

• 270, 13: 4 2 should be 4" 

• 271, 7: one of the (pk)\s should be a (qk)\ 
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